Dynamic buckling and fragmentation in brittle rods 
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We present experiments on the dynamic buckling and fragmentation of slender rods axially im- 
pacted by a projectile. By combining the results of Saint- Venant and elastic beam theory, we derive 
a preferred wavelength A for the buckling instability, and experimentally verify the resulting scal- 
ing law for a range of materials including teflon, dry pasta, glass, and steel. For brittle materials, 
buckling leads to the fragmentation of the rod. Measured fragment length distributions show two 
clear peaks near A/2 and A/4. The non-monotonic nature of the distributions reflect the influence 
of the deterministic buckling process on the more random fragmentation processes. 



Long, thin supports are ubiquitous in natural and en- 
gineered load bearing structures, from spider legs to the 
steel struts of a skyscraper 0,0- A single rod will buckle 
if too much force is applied along its axis, which can lead 
to the catastrophic failure of the structure. The buckling 
instability is seen at all sizes, from pole vaulting to 
protein microtubules confined in vesicles [4 and carbon 
nanotube atomic force microscope probes 0. While the 
classic Euler buckling of a rod is due to a static axial load 
a different physical process occurs when the stress is 
applied suddenly, as during impact 0, @, . I n this Let- 
ter, we show that this "dynamic buckling" obeys a 
simple scaling law, which we derive by combining the ap- 
proach of Saint- Venant with the classical theory of elas- 
tic rods For brittle rods, buckling often leads to 
breaking, for which we find the distribution of fragment 
lengths displays a unique non-monotonic shape reflecting 
the primary buckling instability. 

The process of dynamic buckling and subsequent frag- 
mentation is illustrated in Fig.^ in which a falling weight 
strikes an upright brittle rod (dry pasta) . Within a frac- 
tion of a millisecond after impact, a sinusoidal pertur- 
bation appears (Fig. ^Jj), much different than the half- 
wavelength seen in Euler buckling. A few tenths of a 
millisecond later, the pasta has buckled appreciably and 
begins to shatter (Fig.^i). This imparts angular momen- 
tum of alternating signs to the fragments, which rotate 
and scatter (Fig. ^l-f).Our experimental setup consists 
of a simple metal holder for the rod, and a pneumatic 
cannon in which a pressure reservoir at 60 psi delivers an 
impulse to a steel cylindrical projectile (1.46 cm diam- 
eter, 24.9 or 10.0 g) held at the end of a 1.5 m acrylic 
tube with two magnets. The holder rests on a steel plate 
in a sandbox which serves as a shock absorber to stop 
the projectile. The buckling and fragmentation dynamics 
were imaged by a high speed digital video camera (Phan- 
tom v5.0), capable of capturing up to 62,000 frames per 
second. The speed of the projectile was also measured 
just before impact using the video system. The materi- 
als used included dry pasta (p — 1.5 g/cm 3 , d = 1.1 mm 




FIG. 1: The dynamic buckling and fragmentation of dry 
pasta (d = 1.9 mm, L — 24 cm) just after the impact of 
an aluminum projectile at velocity Uo = 3.5 m/s (interval 
between pictures: 236 /is). 



and 1.9 mm) ^]|, borosilicate glass (p = 2.4 g/cm 3 , d = 
2.0 mm), type 303 stainless steel (p = 7.9 g/cm 3 , d = 
1.6 mm), and teflon (PTFE) (p = 2.2 g/cm 3 , d = 2.0 
mm), all with circular cross-sections and lengths ranging 
from 14-29 cm. We also used strips of plastic and pa- 
per, which had rectangular cross-sections. The Young's 
modulus (E) was measured for the pasta by a clamped 
beam resonance technique , averaging the fundamental 
frequency for five different pieces; we find E ~ 2.9 GPa. 
Standard tabulated values of E were used for the other 
materials: 62 GPa (glass), 200 GPa (stainless steel), and 
0.5 GPa (PTFE). 

The first manifestation of dynamic buckling is the 
onset of undulations with a well-defined wavelength A. 
Our initial experiments involved simply dropping a metal 
weight from a given height, so that the impact velocity 
was fixed; the material, length, and thickness of the rod 
varied. We find that A is strongly dependent on the small- 
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est dimension of the rod (its thickness or diameter d), 
and apparently independent of the length L; in contrast 
for Euler buckling, A = 2L. For a fixed impact veloc- 
ity Uo = 3.5 m/s, measurement of the initial sinusoidal 
perturbation for a wide variety of materials (including 
paper, plastic, and glass) indicates that A ~ d as shown 
in Fig. seemingly independent of the elastic proper- 
ties of each material. A more systematic study, however, 
shows that the material properties do indeed play a role. 
By varying the impact speed Uq from 0.7 to 30 m/s us- 
ing the pneumatic cannon, and widening our study to 
include materials with high and low sound speeds (stain- 
less steel, c = 5020 m/s, and teflon, c = 470 m/s), we 

— 1/2 

find that A ~ U , with a prefactor which varies with 
the material (see Fig. |3^>) ■ 

The sequence of events leading to the buckling of the 
rod begins just after impact, as a stress wave penetrates 
into the material at the speed of sound. In the com- 
pressed region, the rod becomes elastically unstable and 
buckles due to the axial stress. We describe the destabi- 
lization by recalling the dynamical equation for the lat- 
eral dis plac ement £(x,i) of a thin straight elastic beam 
Ullll SHU under an applied longitudinal force F(x,t) 

where A is the cross-sectional area, E is the Young's 
modulus, and / is the moment of inertia of the area 
(I = 7r<i 4 /64 for a circular cross-section). Standard nor- 
mal mode stability analysis of Eq. Q with a constant 
force F(x, t) = Fq yields the fastest growing unstable 
mode, with a wavelength A and associated growth rate 



FIG. 3: The measured dynamic buckling wavelength for var- 
ious materials versus (a) smallest dimension and (b) impact 
speed. The line in (a) corresponds to the scaling A ~ d and 
in (b) corresponds to the scaling A ~ I7 ~ 1/ ' 2 . 



r b 1 given by 

A 7T 3 / 2 [Effi _ x 8F 

where c = y/E/p is the speed of sound. 

The instability is driven by the applied force (Fq), 
which is the elastic response of the material to impact. 
In order to proceed we need to obtain this force from 
the initial momentum of the projectile and the proper- 
ties of the rod. The problem of stress in a rod fixed at 
one end (x = L) and struck longitudinally at the other 
(x = 0) by a mass M with initial speed Uq was solved by 
Saint- Vcnant in 1883 6j. In the linear elasticity limit, the 
longitudinal displacement w(x, t) obeys a wave equation 
dttw — c d xx w and the boundary condition at the impact 
point is the continuity of the impulse Mduw = —EAd x w 
at x = |fj]. The resulting expression for the stress pro- 
file along the rod for < x < ct is then 

where L p = Mj (pA) is the length of the rod which has a 
mass equal to the projectile. Eq.Oflis valid until the front 
has reached the fixed end of the rod (x = L) 13]. After 
that, for t > L/c, the wave rebounds, and the compres- 
sion in the rod is the superposition of the incoming and 
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reflected waves. The full solution of the reflected pulse 
problem can be solved with Saint- Venant's "continuing 
equation", in which the solution w(x,t) during the nth 
time range nL/c < t < (n + l)L/c is obtained from the 
solution during the preceding time range |a,|£|. However, 
to a good approximation, the stress doubles after the first 
rebound. 

Of crucial importance to the shape of the compression 
wave is the characteristic length L p . Typically in our 
experiments M ~ 25 g, which for our experiments gives 
L p = 2 — 4m, about an order of magnitude longer than 
the rod itself! Thus the exponential factors in Eq. (jJJ are 
of order unity, and the compression wave is a step pulse 
travelling at speed c, with magnitude 
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(4) 



Inserting this into Eq. @ leads to the observed trends 
shown in Fig. ©: A ~ d\J c/Uq, independent of L. This 
scaling law is valid if the time scale for the buckling Tb < 
L/c (the instability has to develop before the compression 
front has reached the rod end). This occurs for low sound 
speed materials or large L. However, if the pulse does 
reach the opposite end of the rod before the preferred 
mode develops, a reflected pulse is superimposed on the 
initial one, effectively doubling the stress. We therefore 
anticipate that the wavelength and growth rate will be 
given by 



A 7T / c 

d = V2 V 7CV 



and 



(5) 



with 7 = n + 1 for n reflections. Using this scaling law 
to replot the wavelength data for all of the materials, 
diameters and impact speeds, we find excellent quantita- 
tive agreement over two decades (Fig. 0J). We find that 
all data align with curves with either 7 = 1 (no reflec- 
tions) or 7 = 2 (one reflection), and no other adjustable 
parameters. 

High speed video also allows us to track the growth of 
the preferred buckling mode. For dry pasta (d = 1.1 mm, 
L = 22 cm, 7 = 1) with Uq — 3.5 ± 0.1 m/s, the above 
model predicts Tb = 157±7 /xs. The initial increase of the 
buckling amplitude is indeed well fit by an exponential 
with a time constant of 151 ± 8 us, consistent with the 
linear description of the instability development. 

Note finally that Eq. can be understood hcuristi- 
cally as follows: at time t, the compressive wave front 
has travelled a distance ct. The compression of the rod 
is ~ Uat, which implies a compressive strain e ~ Uot/ct = 
Uq/c. From Hooke's law F/A = Ee ("Ut tensio sic vis"), 
we find F ~ EAU Q /c. 

For brittle materials, dynamic buckling leads rapidly 
to the fragmentation of the rod. Because solid fragmen- 
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FIG. 4: The nondimensional buckling wavelength X/d vs the 
normalized impact speed. The 7 factor accounts for the re- 
flection of the initial impulse (see text). The lines correspond 
to the scaling law in Eq. ©, plotted for 7 = 1 and 2. 



tation typically results in a broad distribution of frag- 
ment sizes, with breaks occurring at apparently random 
locations, a statistical approach is usually taken [l4l Il5| . 
The breakup of a rod under impact is known to give rise 
to broad distributions of fragment sizes 0, 0], with 
a highly skewed shape and a long tail, a phenomenol- 
ogy also encountered in other contexts, including liquid 
sprays [l^. However, the precise relation between the 
global statistics of fragmentation and the instabilities 
which produce it remains largely obscure. 

As is evident from Fig. 1, the wavelength A determined 
by the dynamic buckling process leads to a preferred frag- 
ment length of A/2, because the rod is likely to break at 
the points of highest deformation (maximum rod curva- 
ture). This A-dependent breaking process emerges dis- 
tinctly from a more random fragmentation. To investi- 
gate the interplay of these two processes, we studied the 
impact fragmentation of two different kinds of spaghetti, 
using the same projectile with Uq ~ 3.5 m/s (within 3%). 
We shattered 300 pieces of angel hair pasta, where the 
length ranged from 21-22 cm and d = 1.1 mm (within 
2%). After impact, we measured the length of each frag- 
ment larger than about two diameters, chosen in order 
to remain in the one-dimensional fracture regime, lead- 
ing to a distribution of about 1200 fragments. We also 
shattered 200 pieces of spaghetti, with lengths from 22 - 
23 cm and d — 1.9 mm (within 2%). This resulted in a 
distribution of about 1000 fragments. 

The experimental distributions are shown in Fig. [S] 
In both cases two broad peaks are discernable, indicat- 
ing not one but two preferred lengths. These lengths 
are close to A/2 and A/4 (arrows in the figure). High 
speed video images show that the sinusoidal bending of 
the pasta often has a node at the projectile end, from 
which the distance to a maxiumum curvature point would 
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Fragment Lengths. 

FIG. 5: Size distribution normalized by the experimental 
buckling wavelength A for impact fragmentation of pasta at 
U = 3.5 m/s (c ~ 1400 m/s): a) d = 1.1 mm, A = 48 mm), 
and b) d = 1.9 mm, A = 70 mm, vertically offset for clarity. 
The smooth lines are meant to guide the eye. The arrows 
indicate the fragment sizes A/4 and A/2. The inset shows a 
breaking event (captured at At = 76 /is) contributing to both 
the A/2 and A/4 peaks. 



be A/4 (see inset Fig.[3]and also Fig. lc-d). This would 
contribute one A/4 fragment for each impact. The fact 
that both peaks are slightly lower than expected can be 
explained by concentrated breaking events near the ends 
of the main fragments, from which small pieces detach 
with much less rotation than the larger ones (see Fig|TJ. 

Digital video taken at much higher rates (20,000 fps) 
shows that the breaks do not occur simultaneously, but 
they are nearly so (within tens of microseconds, com- 
pared to the instability growth timsecale of about a hun- 
dred microseconds). This suggests that the fragment dis- 
tribution is determined more by stress dynamics before 
breakup than by a sequential, multiplicative breaking of- 
ten invoked in fragmentation models [19|, as first sug- 
gested by Kolmogorov 20] . However the precise rules 
governing the interaction between the primary compres- 
sion wave, the possible rebounds at the rod extremities, 
and the buckled rod shape, remain to be discovered. 

The Euler buckling of a slender rod is one of the sim- 
plest and most general instabilities of a solid, easily ob- 
servable with almost any material. What we have shown 
here is that, using the ideas of Saint- Venant on elastic 
collisions, the dynamic buckling wavelength can be re- 
lated to the impact velocity. Surprisingly, the scaling 
law for the buckling wavelength, which we have demon- 
strated in teflon, pasta, glass, and steel, depends on only 
two properties of the rod: the speed of sound and the 
diameter. Moreover, the fragmentation of a solid under 
impact, usually conceived of as a random process to be 
treated statistically, has been shown here to include the 
imprint of the deterministic buckling process leading to 



breakup. Our simple experiments suggest that, by "pick- 
ing up the pieces", more complicated problems, such as 
magma fragmentation in explosive volcanic eruptions |2lj 
or the crushing patterns observed in carbon nanotubes 
under axial stress j2^ , may still retain information, albeit 
in statistical form, on the initiating physical processes by 
which the impulse of impact is distributed throughout the 
solid before fragmentation. 

We thank R. Geist for experimental assistance and 
J. Duplat for help in understanding Saint- Venant 's so- 
lution. AB acknowledges support from the National Sci- 
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